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A relativistic Green function approach to the inclusive quasielastic (e, e') scattering is presented. 
The single particle Green function is expanded in terms of the eigenfunctions of the nonhermitian 
optical potential. This allows one to treat final state interactions consistently in the inclusive and 
in the exclusive reactions. Numerical results for the response functions and the cross sections for 
different target nuclei and in a wide range of kinematics are presented and discussed in comparison 
with experimental data. 
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I. INTRODUCTION 

The inclusive electron scattering in the quasielastic re- 
gion addresses to the one-body mechanism as a natu- 
ral interpretation. However, when the experimental data 
of the separation of the longitudinal and transverse re- 
sponses were available, it became clear that the expla- 
nation of both responses necessitated a more compli- 
cated frame than the single particle model coupled to 
one-nucleon knockout. 

A review till 1995 of the experimental data and their 
possible explanations is given in Ref. [1]. Thereafter, 
only few experimental papers were published [2, 3]. 
New experiments with high experimental resolution are 
planned at JLab [4] in order to extract the response func- 
tions. 

From the theoretical side, a wide literature was pro- 
duced in order to explain the main problems raised by the 
the separation, i.e., the lack of strength in the longitudi- 
nal response and the excess of strength in the transverse 
one. The more recent papers are mainly concerned with 
the contribution to the inclusive cross section of meson 
exchange currents and isobar excitations [5-7], with the 
effect of correlations [8, 9], and the use of a relativistic 
frame in the calculations [7]. 

At present, however, the experimental data are not 
yet completely understood. A possible solution could 
be the combined effect of two-body currents and tensor 
correlations [8, 10, 11]. 

In this paper we want to discuss the effects of final 
state interactions in a relativistic frame. Final state in- 
teractions are an important ingredient of the inclusive 
electron scattering, since they are essential to explain 
the exclusive one-nucleon emission, which gives the main 
contribution to the inclusive reaction in the quasielastic 
region. The absorption in a given final state due, e.g., to 
the imaginary part of the optical potential, produces a 
loss of flux that is appropriate for the exclusive process, 
but inconsistent for the inclusive one, where all the al- 
lowed final channels contribute and the total flux must 
be conserved. 

This conservation is preserved in the Green function 



approach considered here, where the components of the 
nuclear response are written in terms of the single particle 
optical model Green function. This result was originally 
derived by arguments based on the multiple scattering 
theory [12] and successively by means of the Feshbach 
projection operator formalism [13-16]. The spectral rep- 
resentation of the single particle Green function, based 
on a biorthogonal expansion in terms of the eigenfunc- 
tions of the nonhermitian optical potential, allows one to 
perform explicit calculations and to treat final state inter- 
actions consistently in the inclusive and in the exclusive 
reactions. Important and peculiar effects are given, in 
the inclusive reactions, by the imaginary part of the op- 
tical potential, which is responsible for the redistribution 
of the strength among different channels. 



In a previous paper of ours [15] the approach was used 
in a nonrelativistic frame to perform explicit calcula- 
tions of the longitudinal and transverse inclusive response 
functions. The main goal of this paper is to extend the 
method to a relativistic frame and produce similar re- 
sults. Although some differences and complications are 
due to the Dirac matrix structure, the formalism follows 
the same steps and approximations as those developed in 
the nonrelativistic frame of Refs. [15, 16]. The numerical 
results obtained in the relativistic approach allow us to 
check the relevance of relativistic effects in the kinemat- 
ics already considered in Ref. [15] and can be applied 
to a wider range of kinematics where the nonrelativistic 
calculations are not reliable. 



In Sec. II the hadron tensor of the inclusive process 
is expressed in term of the relativistic Green function, 
which is reduced in Sec. Ill to a single particle expres- 
sion. The problem of antisymmctrization is discussed in 
Sec. IV. In Sec. V the Green function is calculated in 
terms of the spectral representation related to the opti- 
cal potential. In Sec. VI the results of the calculation 
are reported and compared with the experimental data. 
Summary and conclusions are drawn in Sec. VII. 
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II. THE GREEN FUNCTION APPROACH 
A. Definitions and main properties 

In the one photon exchange approximation the inclu- 
sive cross section for the quasielastic (e, e') scattering on 
a nucleus is given by [1] 



a inc = K (2e L i? L + R T ) , 
where if is a kincmatical factor and 



(1) 



(2) 



measures the polarization of the virtual photon. In Eq. 
(2) ■&(, is the scattering angle of the electron, Q 2 = 
q 2 — u! 2 , and — (uj, q) is the four momentum transfer. 
All nuclear structure information is contained in the lon- 
gitudinal and transverse response functions, Rl and Rt, 
defined by 



R^{u,q) 
R T (u,q) 



W&{<J,q) + W™(w,q) , 



(3) 



in terms of the diagonal components of the hadron tensor 

WZ(u;,q) = f]|(* f |^(9)l*o>| 2 
J-f 

x 6(E +cj- E f ) . (4) 

Here is the nuclear charge-current operator which con- 
nects the initial state | ^o) of the nucleus, of energy 
Eq, with the final states | \&f), of energy Ef, both eigen- 
states of the (A + l)-body Hamiltonian H . The sum 
runs over the scattering states corresponding to all of the 
allowed asymptotic configurations and includes possible 
discrete states. As made for | in the following the de- 
generacy indexes will be omitted whenever unnecessary. 
The ground state | ^o) is assumed to be nondegener- 
ate. In order to avoid complications of little interest in 
the present context, we neglect recoil effects and consider 
only point-like nucleons, without distinguishing between 
protons and neutrons. Unless stated otherwise, the wave 
functions are properly antisymmetrized. 

The hadron tensor of Eq. (4) can equivalently be ex- 
pressed as 



-ilm(*o I Jrt(q)G(Et)J»{q) | *„> 



(5) 



where Ef = Eq + uj and G(Ef) is the Green function 
related to H, i.e., 



G(Ef) 



1 



Ef - H + irj 



(6) 



Here and in all the equations involving G the limit for 
rj — > +0 is understood. It must be performed after calcu- 
lating the matrix elements between normalizable states. 



In this paper the interest is focused on relativistic wave 
functions for initial and final states. Therefore, the (A + 
l)-body Hamiltonian H is the sum of one nucleon free 
Dirac Hamiltonians and two nucleons interactions Vjj>, 
i.e., 



A+l 



A+l 



H=J^(cx j -p j +P j M) + l V n' 



(7) 



where the 4x4 Dirac matrices, otj and (3j, act on the 
bispinor variables of the nucleon j. No particular as- 
sumption is made on the 4x4 matrix structure of Vjj> . 

In order to express the hadron tensor in terms of sin- 
gle particle quantities, the same approximations as in the 
nonrelativistic case [15] are required. The first one con- 
sists in retaining only the one-body part of the charge- 
current operator J M . Thus, we set 



A+l 



■n</) = £jf(<?), 



(8) 



i=i 



where j? acts only on the variables of the nucleon i. By 
Eq. (8), one can express the hadron tensor as the sum of 
two terms, i.e., 



WZ(u,q) = W^(LO,q) 



W c Z(w,q) 



(9) 



where W^(cu,q) is the incoherent hadron tensor [17], 
which contains only the diagonal contributions j^Gjtf, 
whereas the coherent hadron tensor Wf*£(u},q) gathers 
the residual terms of interference between different nu- 
cleons. As the incoherent hadron tensor, also W^{uo,q) 
can be expressed in terms of single particle quantities 
(see Sect. 9 of Ref. [16]), but for the transferred mo- 
menta considered in this paper we can take advantage of 
the high-g approximation [18] and retain only W^(u, q). 
This term can be further simplified using the symmetry of 
G for the exchange of nucleons and the antisymmetriza- 
tion of | \&o)- Therefore we write 



A+l 



x Im(* | jrt(q)G(Et)j»(q) | *„) ,(10) 

where j^{q) is the component of J^{q) related to an ar- 
bitrarily selected nucleon. Due to the well-known com- 
pleteness property, i.e., 

J- J AE(GHE)-G(E))=l, (11) 

the incoherent hadron tensor fulfills the energy sum rule 

j <W"(w,g) = (A+l)(V | jrt(q)j"(q) | *„) • (12) 

B. Projection operator formalism 

This formalism yields an expression of the incoherent 
hadron tensor of Eq. (10) in terms of eigenfunctions and 
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Green functions of the optical potentials related to the 
various reaction channels. Apart from complications due 
to the Dirac matrix structure, we follow the same steps 
and approximations as in the nonrelativistic treatment 
[12, 13, 15, 16]. 
Let us decompose H as 



H = a-p + (3M + U + H R 



(13) 



where a-p+f3M is the kinetic energy of an arbitrarily se- 
lected nucleon, U is the interaction between this nucleon 
and the other ones, and Hr is the residual Hamiltonian 
of A interacting nucleons. Such a decomposition cannot 
be performed in the physical space of the totally antisym- 
metrizcd (A + 1) nucleon wave functions. Therefore, we 
must operate in the Hilbert space Tt of the wave func- 
tions which are antisymmetrized only for exchanges of 
the nucleons of H R . This treatment is presented here 
only for sake of simplicity. In Sect. IV we shall discuss 
its physical drawbacks and outline the necessary changes. 

Let | n) and | e) denote the antisymmetrized eigenvec- 
tors of Hr related to the discrete and continuous eigen- 
values e n and e, respectively. We introduce the operators 
P n , projecting onto the n-channel subspace of H 7 and Q n , 
projecting onto the orthogonal complementary subspace, 
i.e., 

P n = ^ / dr | ra;n)(n;ra \ , 

a J 

Qn = l-Pn- (14) 

Here | ra;n) is the unsymmctrized vector obtained from 
the tensor product between the discrete eigenstate | n) 
of Hr, and the orthonormalized eigenvectors | ra) (a = 
1, 2, 3, 4) of the position and the spin of the selected nu- 
cleon. The eigenvectors | ra) have been chosen only for 
sake of dcfinitcness, as every complete orthonormalized 
set of single nucleon vectors would define the same opera- 
tors P n . Apart from minor differences due to the present 
relativistic context, P n and Q n are the projection op- 
erators of the Feshbach unsymmetrized formalism [19]. 
Note, for later use, the relations 



[P n ,a-p + f3M] = 

HrPu = £ n Pn ■ 



(15) 



Moreover, we introduce the projection operator onto the 
continuous channel subspace, i.e., 



dr | ra; e)(e; ra 



(16) 



Due to the completeness of the set {| ra; n), | ra; e)}, one 
has 



(17) 



Then, we insert Eq. (17) into Eq. (10) disregarding the 
contribution of P c . This approximation, which simplifies 



the calculations, is correct for sufficiently high values of 
the transferred momentum q. Thus, the hadron tensor 
of Eq. (10) can be expressed as the sum 

W^(w,q) = W^(w,q) + W^(w,q) , (18) 

of a direct term 

n 

Wr(u,q) = -A±li m (^ I jrt(q)P n G(E { )P n 

x f{q) | *o> , (19) 

and of a term 

W^(u>,q) = (",</) . 

n 

Wr(",q) = -^-lm{* \jrt{q)P n G{E t )Q n 

x f{q) | *o> , (20) 

which gathers the contributions due to the interference 
between the intermediate states | ra; n) related to differ- 
ent channels. 

We note that the interference term does not contribute 
to the energy sum rule. In fact Eq. (11) yields 



/ 



x F(q) | * > = • (21) 



Thus, the full contribution to the sum rule of the inco- 
herent hadron tensor is given only by the direct term, 
i.e., 



/ 



du;W^(LO,q) = J dujW^(LU,q) 

= (A + 1)<* | f\q) ]T P n f{q) | *o> , (22) 



which, as a pure consequence of the omission of the con- 
tinuous channels described by P c , is smaller than the 
value of Eq. (12). 



III. SINGLE PARTICLE EXPRESSION OF THE 
HADRON TENSOR 

A. Single particle Green functions 

For the time being, we disregard the effects of inter- 
ference between different channels and consider only the 
direct contribution to the hadron tensor of Eq. (19). The 
matrix elements of P n G(E)P n in the basis | ra; n) define 
a single particle Green function Q n {E) having a 4 x 4 
matrix structure, i.e., 

(ra | g n (E) | rV) = (n; ra | G(E + e n ) \ r'a';n) . (23) 



4 



Note that here the energy scale is in accordance with Ref. 
[16] and differs from Ref. [15]. 

The self-energy of G n (E) is determined following the 
same steps used by Feshbach to determine the optical 
potential from the Schrodingcr equation [19]. One starts 
from the relation 

{E - a • p - 0M - H R - U + e n + it]) (P n + Q n ) 

xG(E + e n )P n = P n , (24) 

projects both sides by P n and then by Q n , uses Eqs. (15), 
resolves Q n GP n in terms of P n GP n , and finally obtains 



{E - a - p - 0M -V n (E) + ir/) 
xP n G(E + s n )P n = P n , 



(25) 



with 

V n (E) = P n UP n 

+ PnUQr 



1 



E - Q n HQ n + e n + ii] 



>UP n .(26) 



Using Eq. (14) for P n and considering the matrix ele- 
ments in the basis | ra; n) of both sides of Eq. (25), one 
has 



Q n {E) 



1 



(27) 



E-h n (E)+ir, ' 
where 

h n (E) = a-p + 0M + V n (E) , (28) 

and V n {E) has the 4x4 matrix structure defined by 

(ra | V n (E) | rV> = (n; ra | V n (E) | r'a'; n) . (29) 

Thus, h n (E) is the self-energy of the Green function 
Gn(E) and V n (E) is the related mean field. Using the 
same arguments as in the nonrelativistic case, one finds 
that V n (E) is the unsymmetrized Feshbach optical po- 
tential [19], related to the channel n, for the relativistic 
Hamiltonian H . 

Using the first Eq. (14)andEq. (23), the direct hadron 
tensor W^(lu, q) of Eq. (19) becomes 



1 



x f{q) | Vn) , 



(30) 



where the initial state | <p n ), normalized to 1, is repre- 
sented by the bispinor defined by the matrix elements 



(ra | ip n ) 



A+l 



(n;ra | *„} 



A„ is the related spectral strength [20] 

\ n = (A+l)J2 f dr| (n;ra|* > | 2 



(31) 



(32) 



with 

]T A n ~ A + 1 , (33) 

n 

and the symbol (/ | g) denotes the scalar product 

(f\g)=J2[ drf(ra)g(ra) . (34) 

In Eq. (30) the hadron tensor is expressed in terms of 
single particles quantities. As in the nonrelativistic case, 
| ip n ) are the eigenstates of the optical potential, i.e., 



(a-p + pM + V n (E Q -£„)) | ip n ) 

= (E - £ n ) | ip n ) . 



(35) 



If \^e) is the eigenstate of H corresponding asymptoti- 
cally to a nucleon, of momentum k, colliding with a tar- 
get nucleus in the bound state | e„), the single particle 
vectors | Xn(E — £„)) representing the elastic scattering 
wave functions (n; ra | ^e) are eigenstates of the same 
optical potential, i.e., 

(a ■ p + 0M + V n (E - £„)) | Xn (E - £„)) 

= (E-e n ) | X n(E-e n )) . (36) 



Since E is the total energy \/fc 2 + M 2 + e n , the argument 
of V„ is the kinetic energy (including the rest mass) of 
the emitted nucleon. 



B. Interference hadron tensor 

The problem of expressing the interference hadron ten- 
sor in a one-body form is treated in Ref. [15] in 
the nonrelativistic context. It is argued that the contri- 
bution of can be included into the direct hadron 
tensor by the simple replacement 



g n (E)^gf(E) y/i-v' n (E)g n (E) 



x Jl-V n (E) 



(37) 



where V n (E) is the energy derivative of the Feshbach 
optical potential. 

In Ref. [21] the problem is considered anew from a 
rigorous point of view. The interference hadron tensor is 
expressed exactly as a series involving energy derivatives 
of V„ (E) , of increasing order, plus a residual term which 
cannot be reduced to a single particle form. The series is 
expected to fastly converge near the quasielastic peak and 
at intermediate energies. It is argued that in this region 
of momenta and energies the residual term is negligible. 
Thus, one recovers the result of Eq. (37) and second 
order corrections which do not seem to give a sizable 
contribution. 

Neither the treatment nor the conclusions change if 
one considers the relativistic Hamiltonian H . Thus, for 
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the hadron tensor of Eq. (18) we use the approximated 
expression obtained from Eq. (30) with the replacement 
(37): 

7T * — ' 

n 

x f(q) | <p n ) . (38) 

Since the interference term of Eq. (20) has no in- 

fluence on the energy sum rule of the total hadron tensor 
of Eq. (18), a natural question arises whether the ap- 
proximation leading to Eq. (38) may change the sum 
rule. Actually, one can observe that in Eq. (37) Gn(E) is 
modified by factors which change neither its properties 
of analyticity in the energy complex plane nor its high 
energy behavior. This fact is used in Rcf. [15] to prove 
the relation 

-^J dE Im (™ I Sf(E) | rV) 
= ~ J dE Im(ra | Q n (E) | r'a') 
= S (r - r') 5 aa < • (39) 

Therefore, the energy sum rule obtained from Eq. (38) 
is exactly the same as in Eq. (22), i.e., the correct sum 
rule of the incoherent hadron tensor, apart from the con- 
tribution of the continuous channels. 



C. Excited states of the residual nucleus 

As neither microscopic nor empirical calculations are 
available for the optical potential V„ associated with the 
excited states | e n ), a common practice relates them to 
the ground state potential Vo by means of an appropriate 
energy shift. Here, as in Ref. [15], we use the kinetic en- 
ergy prescription for the shifts (see Sect. 5 of Ref. [16]), 
naturally suggested by the plane wave impulse approx- 
imation. Such a prescription preserves the value of the 
kinetic energy (including the rest mass), directly related 
to the value of the optical potential variable in the energy 
scale adopted here. Therefore we set 

V n (E) ~ V (E) , (40) 

which implies 

g n (E) ~ g (E) . (4i) 

Using these approximations in Eq. (38), we write 
W»»(u>,q) = --Y j \ n \m{y n \jrt{q)Qf{E i -e n ) 

7T * — ' 

n 

x j"(q) I <p n ) ■ (42) 

These approximations do not change the energy sum rule 
of W^{u>,q). 



IV. ANTISYMMETRIZATION 

For sake of simplicity the treatment of Sects. II and 
III is based on the unsymmetrized projection operator 
P n defined in Eq. (14), leading to the Green function 
Q n of Eq. (23). In this Section we examine the draw- 
backs of this formulation and the possible alternatives. 
On the mathematical ground, Q n deserves the name of 
Green function since it fulfills the sum rule (39), which is 
a qualifying property. Moreover, and intimately related, 
Q n is an invertible operator on the whole Hilbert space 
L 2 (M 3 ). Hence, its self-energy is not affected by any un- 
due restriction of domain and by the related mathemat- 
ical troubles. 

Notwithstanding, the optical potential V n related to 
Q n suffers from the drawback of having spurious eigen- 
functions. In fact H has both antisymmetrized and un- 
symmetrized eigenvectors | ^e) and the latter ones gen- 
erate eigenfunctions (n; ra | ^>e) of V n (E — e n ) which 
have no physical meaning. No tool exists to make a dis- 
tinction inside this unphysical degeneracy. Besides, it is 
apparent that V n cannot be compared with any empirical 
optical model potential. 

The remedy is a treatment based on projection op- 
erators onto antisymmetrized states, although their in- 
clusion into the hadron tensor is more laborious. Two 
approaches are available. 

a) The first one (see Subsect. 3.2 of Rcf. [16] and 
Appendix B of Ref. [15]) uses the Feshbach projection 
operator onto antisymmetrized states, i.e., 

P n = E / drdr'atjn) 

a. a' 

x (Ta\(l-K n )- 1 \T f a')(n\a T , a , , (43) 

where a\ a creates a nucleon in the state | ra) and K n is 
the one-body density matrix defined by 

(ra | K n | r'a') = (n \ al a ,a ra \ n) . (44) 

The results of the previous Section remain true with the 
replacement 

(ra | Qn(E) \ rV) -> (ra | Q F n (E) \ rV) 

x (s6 | (1 - Kn)' 1 I r'a') , (45) 

where Q F is the Green function related to the "sym- 
metrized" Feshbach optical potential [22]. In this 
approach the spurious degeneracy disappears, since one 
operates in a Hilbert space of antisymmetrized states, 
but at the price of new drawbacks [16, 23], namely: (1) 
Qn is not fully invertible and so in some cases it gives 
rise to incorrect Dyson equations; (2) both Q F and 
are not symmetric for exchange ra <-> r'a' and therefore 
Vj, is nonhcrmitian below the threshold of the inelastic 
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processes; (3) the usual nonlocal models of potential are 
probably inadequate for V„, which shows a complicate 
nonlocal structure. 

In short, the approach based on disguises nontrivial 
mathematical problems and it is not really useful, since 
this potential bears no close relation with the empirical 
optical model potential. 

b) The second approach, where the above drawbacks 
disappear, is the one of Ref. [16]. It is based on the 
extended projection operator of Ref. [24] 

Pi P+h) = E / dr(at a + a ra ) |n) 

x (« | (4a + O , (46) 

which leads to 

(ra\Q n (E) \v'a')^{va\g^ +h \E) | rV) 
( U I a ™~B it , I 7Z a l'a' I n ) 



+ (« I <£> a > 



E — H + e n + it] 
1 



E + H - e n - irj 



a ra | n) . 



(47) 



g^ +h \E) is the full Green function, including particle 
and hole contributions. It fulfills the sum rule of Eq. 
(39), is fully invertiblc and produces mathematically cor- 
rect Dyson equations. The related mean field V^ +h) (E) 
has no spurious eigenfunctions corresponding to unsym- 
mctrized states and its properties of nonlocality and sym- 
metry make it more easily comparable with the empirical 
optical model potentials. Therefore, we understand that 
in the following equations Q n will denote the full Green 
function Q^ +h ^ of the rclativistic Hamiltonian H. The 
associated mean field V„ is nonlocal as in the nonrel- 
ativistic case and does not conserve the primarily 4x4 
matrix structure of Vjj' . 



V. SPECTRAL REPRESENTATION OF THE 
HADRON TENSOR 



In this Section we consider the spectral representation 
of the single particle Green function which allows prac- 
tical calculations of the hadron tensor of Eq. (42). In 
expanded form, it reads 

W*{w,q) = --yXlmfcn I j»Hq)Jl-V'(E) 

TT L — ' v 



x Q{E)^\-V\E)f{q)\^ n ) , (48) 

where E = Ef — e n . Here and below, the lower index 
is omitted in the Green functions and in the related 
quantities. According to the discussion of the previous 
Section, we understand that Q is the full particle-hole 
Green function of Eq. (47) and that V is the mean field 



potential related to Q by the equations 
S(E) 



(49) 
(50) 



E - h(E) + irj ' 
h{E) = a p + (3M + V(E) . 

The use of this Green function does not change the ex- 
pressions of the normalized initial states | ip n ) and of the 
related spectroscopic factors A„, defined in Eqs. (31) and 
(32), respectively. Equivalently, they can be written as 



(ra 



An 



dr 



a | *0> , 

(n | a ra 



(51) 
(52) 



Due to the complex nature of V{E) the spectral represen- 
tation of Q (E) involves a biorthogonal expansion in terms 
of the eigenfunctions of h(E) and h){E). We consider 
the incoming wave scattering solutions of the eigenvalue 
equations, i.e., 



h\E)) 
-HE)) 



(53) 
(54) 



The choice of incoming wave solutions is not strictly nec- 
essary, but it is convenient in order to have a closer com- 
parison with the treatment of the exclusive reactions, 
where the final states fulfill this asymptotic condition and 
are the eigenfunctions | \e °f h^(E). 

The eigenfunctions of Eqs. (53) and (54) satisfy the 
biorthogonality condition 



(x < f ) (E)\x ( ->(E))= 6 (£-£') , 

and, in absence of bound eigenstates, the completeness 
relation 



(55) 



d£ 



M 



(e))(x { £ -\e) 



(56) 



where the nucleon mass M is the threshold of the con- 
tinuum of h(E). 

Eqs. (55) and (56) are the mathematical basis for the 
biorthogonal expansions. The contribution of possible 
bound state solutions of Eqs. (53) and (54) can be disre- 
garded in Eq. (56) since their effect on the hadron tensor 
is negligible at the energy and momentum transfers con- 
sidered in this paper. 

Inserting Eq. (56) into Eq. (49) and using Eq. (54), 
one obtains the spectral representation 



G{E) 



1 



(E)) 
(X £ '\E) 



E-£ + ir] 
Therefore, Eq. (48) can be written as 
1 x - 



(57) 



W^(lo,q) 



Im 



d£ 



M 



E { - 



£ + ir] 



x TW{£,Et-e n ) 



(58) 
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where 



Tr(£,E) = \ n {v n \f\q)\ll-V{E)\x { -\E)) 

x <x<f\E) | y/l -V(E)j»(q) \ Vn ) . (59) 

The limit for r\ — ► +0, understood before the integral 
of Eq. (58), can be calculated exploiting the standard 
symbolic relation 



1 



lim 

r;^0 E - £ +17] 



1 



-Mr£(.E-£) , (60) 



where P denotes the principal value of the integral. 
Therefore, Eq. (58) reads 



ReT^{Ei-s n ,Ei-s n ) 



k Ja 



d£- - 

Ef — e n — £ 



lmTW{£,E f -e n ) 



(61) 



which separately involves the real and imaginary parts of 

^ n • 

Some remarks on Eqs. (59) and (61) are in order. Let 
us examine the expression of T^(£,E) at £ = E = 
Ei — e n for a fixed n. This is the most important case 
since it appears in the first term in the right hand side 
of Eq. (61), which gives the main contribution. Disre- 
garding the square root correction, due to interference 
effects, one observes that in Eq. (59) the second matrix 
element (with the inclusion of V^n) is the transition am- 
plitude for the single nucleon knockout from a nucleus in 
the state | ^o) leaving the residual nucleus in the state 
| n). The attenuation of its strength, mathematically 
due to the imaginary part of V , is related to the flux 
lost towards the channels different from n. In the inclu- 
sive response this attenuation must be compensated by 
a corresponding gain due to the flux lost, towards the 
channel n, by the other final states asymptotically origi- 
nated by the channels different from n. In the description 
provided by the spectral representation of Eq. (61), the 
compensation is performed by the first matrix element 
in the right hand side of Eq. (59), where the imaginary 
part of V has the effect of increasing the strength. Similar 
considerations can be made, on the purely mathematical 
ground, for the integral of Eq. (61), where the ampli- 
tudes involved in have no evident physical meaning 
as £ ^ Ef — e„. We want to stress here that in the 
Green function approach it is just the imaginary part of 
V which accounts for the redistribution of the strength 
among different channels. 

The matrix elements in Eq. (59) contain the mean field 
V(E) and its hermitian conjugate V^(-E), which are non- 
local operator with a possibly complicated matrix struc- 
ture. Neither microscopic nor empirical calculations of 
V(E) are available. In contrast, phcnomcnological opti- 
cal potentials are available. They are obtained from fits 



to experimental data, are local and involve scalar and 
vector components only. The necessary replacement of 
the mean field by the empirical optical model potential 
is, however, a delicate step. 

In the nonrelativistic treatment of Rcfs. [15, 21] this 
replacement is justified on the basis of the approximated 
equation (holding for every state | ip)) 



Im(V | ^l-V{E)g(E)^l-V{E) | V) 
; Im<V | Jl~V' L (E)gL(E)Jl-V L (E) | V) , (62) 



where Vi,{E) is the local phase-equivalent potential iden- 
tified with the phcnomcnological optical model potential 
and Gl(E) is the related Green function. In Rcf. [21] the 
proof of Eq. (62) is based on two reasons: (1) a model 
of V{E) commonly used in dispersion relation analyses; 
(2) the combined effect of the factor \J\ — V'(E) and 
of the Perey factor, which connects the eigenfunctions 
of V(E) and Vl(E). We stress that it is just the fac- 
tor \J\ — V'(E), introduced to account for interference 
effects, which allows the replacement of V(E) by Vl(-E). 

Although the Perey effect is not sufficiently known for 
the Dirac equation, we have a reasonable confidence that 
Eq. (62) holds also in the present relativistic context. 
Therefore, we insert Eq. (62) into Eq. (48). Then, all 
the developments of this Section can be repeated with 
the simple replacement of V(E) by Vl(E). 



VI. RESULTS AND DISCUSSION 

The cross sections and the response functions of the 
inclusive quasielastic electron scattering are calculated 
from the single particle expression of the coherent hadron 
tensor in Eq. (61). After the replacement of the 
mean field V(E) by the empirical optical model potential 
Vl (E) , the matrix elements of the nuclear current oper- 
ator in Eq. (59), which represent the main ingredients 
of the calculation, are of the same kind as those giv- 
ing the transition amplitudes of the electron induced nu- 
cleon knockout reaction in the relativistic distorted wave 
impulse approximation (RDWIA) [25]. Thus, the same 
treatment can be used which was successfully applied to 
describe exclusive (e, e'p) and (j,p) data [25, 26]. 

The final wave function is written in terms of its upper 
component following the direct Pauli reduction scheme, 
i.e., 



x { £ -\e)= 



1 

M + £ + & (E) - Vl (E) 



f+ 



(63) 



where S(E) and V(E) are the scalar and vector energy- 
dependent components of the relativistic optical poten- 
tial for a nucleon with energy E [27]. The upper com- 
ponent, 'Sf+j is related to a Schrodingcr equivalent wave 
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function, <j>f, by the Darwin factor, i.e., 
* f+ = \lD\{E)* t , 
D £ (E) = 1 



S(E) - V{E) 



M + £ 



(64) 
(65) 



<&f is a two-component wave function which is solution of 
a Schrodinger equation containing equivalent central and 
spin-orbit potentials obtained from the scalar and vector 
potentials [28, 29]. 

As no relativistic optical potentials are available for 
the bound states, then the wave function cp n is taken 
as the Dirac-Hartree solution of a relativistic Lagrangian 
containing scalar and vector potentials [30, 31]. 

Concerning the nuclear current operator, no unam- 
biguous approach exists for dealing with off-shell nuclc- 
ons. In the present calculations we use the cc2 expression 
of the 1-body current [32-34] 



Jcc2 



WV + i^MQVV, (66) 

where k is the anomalous part of the magnetic moment, 
Fi and F 2 are the Dirac and Pauli nucleon form factors, 
which are taken from Ref. [35], and a» v = (i/2) [7^,7"]. 

Current conservation is restored by replacing the lon- 
gitudinal current by [33] 



J L = J z = 



J 



(67) 



The calculations have been performed with the same 
bound state wave functions and optical potentials as in 
Rcfs. [25, 26], where the RDWIA was able to reproduce 
(e, e'p) and (7,2?) data. 

The relativistic bound state wave functions have been 
obtained from the code of Rcf. [30], where relativistic 
Hartree-Bogoliubov equations are solved in the context 
of a relativistic mean field theory to reproduce single- 
particle properties of several spherical and deformed nu- 
clei [31]. The scattering state is calculated by means of 
the energy- and mass number-dependent EDAD1 com- 
plex phenomenological optical potential of Rcf. [27], 
which is fitted to proton elastic scattering data on several 
nuclei in an energy range up to 1040 MeV. 

In the calculations the residual nucleus states | n) are 
restricted to be single particle one-hole states in the tar- 
get. A pure shell model is assumed for the nuclear struc- 
ture, i.e., we take a unitary spectral strength for each 
single particle state and the sum runs over all the occu- 
pied states. 

The results presented in the following contain the con- 
tributions of both terms in Eq. (61). The calculation of 
the second term, which requires integration over all the 
eigenfunctions of the continuum spectrum of the optical 
potential, is a rather complicate task. This term was ne- 
glected in the nonrelativistic investigation of Ref. [15], 
where its contribution was estimated to be very small. 
In the present relativistic calculations the effect of this 
term can be significant and it is therefore included in the 
results. 



C(e,e') 




250 300 

to [MeV] 



„ 0.03 




250 300 

oj [MeV] 

FIG. 1: Longitudinal (upper panel) and transverse (lower 
panel) response functions for the 12 C(e, e') reaction at q = 400 
MeV/c. Solid and dotted lines represent the NLSH results 
with and without the inclusion of the factor in Eq. (68), re- 
spectively. Dashed lines give the result without the integral in 
Eq. (61). Dot-dashed lines are the contribution of integrated 
single nucleon knockout only. The data are from Ref. [36] . 

A. The i?L and Rt response functions 

The longitudinal and transverse response functions for 
12 C at q = 400 MeV/c arc displayed in Fig. 1 and com- 
pared with the Saclay data [36] . The low energy transfer 
values are not given because the relativistic optical po- 
tentials are not available at low energies. 

The agreement with the data is generally satisfactory 
for the longitudinal response. In contrast, the transverse 
response is underestimated. This is a systematic result of 
the calculations. It may be attributed to physical effects 
which are not considered in the present approach, e.g., 
meson exchange currents. 

The effect of the integral in Eq. (61) is also displayed. 
At variance with the nonrelativistic result, here this con- 
tribution is important and essential to reproduce the ex- 
perimental longitudinal response. 

As explained in Sect. Ill, the contribution arising from 
interference between different channels, see Eqs. (37) and 
(62), gives rise to the factor 



sj 1 - V L (E) = y/l- WW) - V'{E) 



(68) 



We see, however, that here it gives only a slight contribu- 
tion, due to a compensation between the energy deriva- 
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FIG. 2: Longitudinal (upper panel) and transverse (lower 
panel) response functions for the 12 C(e, e') reaction at q — 400 
MeV/c. Solid lines represent the NLSH results, dashed lines 
the NL3 results. Data as in Fig. 1. 



tives S'(E) and V'(E), while in the nonrelativistic cal- 
culation an overall reduction was observed, which was 
necessary to reproduce the data [15]. 

The contribution from all the integrated single nucleon 
knockout channels is also drawn in Fig. 1. It is signif- 
icantly smaller than the complete calculation. The re- 
duction, which is larger at lower values of u, gives an 
indication of the relevance of inelastic channels. 

For the calculations in Fig. 1 the Hartree-Bogoliubov 
equations for the single particle bound states have been 
solved using NLSH choices for the parameters of the rel- 
ativists mean field theory Lagrangian. In Fig. 2 a com- 
parison is shown between the results obtained with this 
choice of parameters and a different choice, i.e., NL3. 
The shapes of the responses calculated with the differ- 
ent bound states show small differences. Their integrals 
must be unchanged, according to the fact that the sum 
rule has to be preserved. 

The longitudinal and transverse response functions for 
12 C at q = 500 and q = 600 MeV/c are displayed in Figs. 
3 and 4, respectively, and compared with the Saclay data 
[36]. The bound state wave function have been obtained 
with the NLSH parametrization. As already found in Fig. 
1 at q = 400 MeV/c, a good agreement with the data 
is obtained in both cases for the longitudinal response, 
while the transverse response is underestimated. Also in 
Figs. 3 and 4 only a slight effect is given by the factor in 
Eq. (68) arising from the interference between different 
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FIG. 3: The same as in Fig. 1, but for q = 500 MeV/c. The 
data are from Ref. [36] . 



channels. The role of the integral in Eq. (61) decreases 
increasing the momentum transfer. At q — 500 MeV/c 
its contribution is smaller than at q = 400 MeV/c, but 
still important to reproduce the experimental longitudi- 
nal response, while at q = 600 MeV/c the effect is neg- 
ligible and the two curves with and without the integral 
overlap. The effect of the inelastic channels, indicated 
in the figures by the difference between the complete re- 
sults and the contribution from all the integrated single 
nucleon knockout channels, is always visible and even siz- 
able, but it decreases increasing the momentum transfer. 

The response functions for 40 Ca at q = 450 MeV/c 
are shown in Fig. 5 and compared with the Saclay [37] 
and the MIT-Bates [3] data. The results obtained with 
the NLSH set of parameters have been plotted, since the 
results with other sets are almost equivalent. The calcu- 
lated response functions are of the same order of mag- 
nitude as the MIT-Bates data, while for the Saclay data 
the longitudinal response is overestimated and the trans- 
verse response underestimated. The factor in Eq. (68) 
produces and enhancement which is minimal but visible 
in the figure. 



B. The inclusive cross section 

Investigation of inclusive electron scattering in the re- 
gion of large q is of great interest to provide informa- 
tion on the nuclear wave functions and excitation and 
decay of nucleon resonances. Several experiments have 
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FIG. 4: The same as in Fig. 1, but for q = 600 MeV/c. The 
data are from Ref. [36]. 



been carried out to explore this region. The separa- 
tion of the longitudinal and transverse components of 
the nuclear response would be very interesting, but it is 
very difficult to perform because of the decreasing of the 
longitudinal-transverse ratio with increasing q. Precise 
measurements over a kinematical range that would al- 
low longitudinal-transverse separation for several nuclei 
are however planned in the future at JLab, where the E- 
01-016 approved experiment [4] will make a precise mea- 
surement in the momentum transfer range 0.55 < q < 1.0 
GeV/c in order to extract the response functions. 

In this Subsection we focus our attention on experi- 
mental cross sections with ui < 300 MeV, since our model 
does not include meson exchange currents and isobar ex- 
citation contributions. 

The calculated inclusive 12 C(e,e') cross section is dis- 
played in Fig. 6 in comparison with the SLAC data [38] 
in a kinematics with a beam energy E e = 2020 MeV 
and a scattering angle of ~ 15°. The bound state wave 
function has been obtained with the NLSH set. A visi- 
ble enhancement is produced by the factor in Eq. (68). 
The effect of the integral in Eq. (61) gives a significant 
reduction which underestimates the data. As in the case 
of the transverse response of Figs. 1-5, the discrepancy 
might be due to two-body mechanisms which are here 
neglected. 

In order to extend our analysis to different kinematics 
and target nuclei, we consider in Fig. 7 the 16 0(e,e') 
inclusive cross section data taken at ADONE-Frascati 
[2] with beam energy ranging from 700 to 1200 MeV and 




250 300 

oj [MeV] 

FIG. 5: Longitudinal (upper panel) and transverse (lower 
panel) response functions for the 40 Ca(e,e') reaction at q = 
450 MeV/c. The Saclay data (open circles) are from Ref. 
[37], the MIT-Bates (black circles) are from Ref. [3]. Line 
convention as in Fig. 1. 



a scattering angle of ~ 32°. The NLSH wave functions 
have been used in the calculations. The agreement with 
data is good in all the considered situations. The integral 
in Eq. (61) produces a reduction which is now essential 
to reproduce the data at 700 MeV, which correspond to 
a momentum transfer of < 400 MeV/c. Its contribution 
can be neglected in the other kinematics, where q ~ 600 
MeV/c. The effect of the factor in Eq. (68) is very small. 



VII. SUMMARY AND CONCLUSIONS 

A relativistic approach to inclusive electron scattering 
in the quasielastic region has been presented. This work 
can be considered as an extension of the nonrelativistic 
many-body approach of Ref. [15]. The components of 
the hadron tensor are written in terms of Green functions 
of the optical potentials related to the various reaction 
channels. The projection operator formalism is used to 
derive this result. An explicit calculation of the single 
particle Green function can be avoided by means of its 
spectral representation, based on a biorthogonal expan- 
sion in terms of the eigenfunctions of the nonhcrmitian 
optical potential V(E) and of its hermitian conjugate. 
The interference betwee n different channels is taken into 
account by the factor \/l — V'(E), which also allows the 
replacement of the mean field V(E) by the phenomeno- 
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FIG. 6: The cross section for the inclusive 12 C(e,e') reaction 
at tf e = 15.02° and E e = 2020 MeV. The data are from SLAC 
[38]. Line convention as in Fig. 1. 
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FIG. 7: The cross section for the inclusive 16 0(e,e') reaction 
at tf e = 32° and E e = 700, 1080, and 1200 MeV. The data 
are from ADONE-Frascati [2]. Line convention as in Fig. 1. 



logical optical potential Vl(E). After this replacement, 
the nuclear response functions are expressed in terms of 
matrix elements similar to the ones which appear in the 
exclusive one nucleon knockout reactions, and the same 
RDWIA treatment [25] can be applied to the calculation 
of the inclusive electron scattering. 

The effects of final state interactions are thus described 
consistently in exclusive and inclusive processes. Both 
the real and imaginary parts of the optical potential must 
be included. In the exclusive reaction the imaginary part 
accounts for the flux lost towards other final states. In 
the inclusive reaction, where all the final states are in- 
cluded, the imaginary part accounts for the redistribu- 
tion of the strength among the different channels. 

All the final states contributing to the inclusive reac- 
tion arc contained in the Green function, and not only 
one nucleon emission. Our calculations for the inclusive 
electron scattering are different from the contribution of 
integrated single nucleon knockout only. The difference 
between the two results is originated by the imaginary 
part of the optical potential. 

The transition matrix elements are calculated using the 
bound state wave functions obtained in the framework of 
a relativistic mean field theory. The direct Pauli reduc- 
tion method is applied to the scattering wave functions. 

Numerical results for the longitudinal and transverse 
response functions of 12 C and 40 Ca have been presented 
in comparison with data in a momentum transfer range 
between 400 and 600 MeV/c. 



The role and relevance of the various effects of final 
state interactions can be different in the relativistic and 
nonrclativistic calculations. This is a consequence of the 
different features of the optical potentials in the two ap- 
proaches. The final effect is however similar and pro- 
duces qualitatively similar results in comparison with 
data. The relativistic frame has anyhow the advantage 
that it can more reliably be applied to a wider range of 
situations and kinematics. 

Our relativistic results confirm that the effects of final 
state interactions are large and essential to reproduce the 
data. The term with the integral, entering the definition 
of the hadron tensor W^^luj.q) in Eq. (61), gives a sig- 
nificant contribution, which is important to improve the 
agreement with data. This result is different from the 
one obtained in the nonrelativistic analysis [15], where 
this term gave only a small contribution and was thus 
neglected in the calculations. We stress that this term is 
due to the imaginary part of the optical potential, which 
thus produces different but important effects in the rela- 
tivistic and nonrelativistic approaches. 

The effects of the integral in Eq. (61) as well as the dif- 
ference between the complete result and the contribution 
of integrated single nucleon knockout, which are both en- 
tirely due to the imaginary part of the optical potential, 
tend to decrease with increasing momentum transfer. 

The factor \J\ — V'(E) is conceptually very important. 
It accounts for interference effects and allows the replace- 
ment of V(E) by Vl(E). In the nonrelativistic analysis 
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of Ref. [15] this factor produced an overall reduction of 
the calculated strength which significantly improved the 
agreement with the experimental longitudinal response 
function. Only a small contribution is given by this fac- 
tor in the present relativistic approach. It generally pro- 
duces a small enhancement of the calculated responses 
that does not change significantly the comparison with 
data. 

Final state interactions have a similar effect on the 
longitudinal and transverse components of the nuclear 
response. In comparison with data, the longitudinal re- 
sponse is usually well reproduced, while the transverse 



response is underestimated. This seems to indicate that 
more complicated effects, e.g., two-body meson exchange 
currents, have to be added to the present single particle 
approach. 

The inclusive cross section for 12 C and 16 O has been 
calculated for momentum transfer < 600 MeV/c. The 
results for 12 C are in agreement with those obtained for 
the response functions. The lack of strength in the de- 
termination of the transverse response results in an un- 
derestimation of the data. A satisfactory agreement is 
obtained for the 16 0(e, e') results. 
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